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1. Introduction 

After Weinberg suggested the application of chiral perturbation theory, a low energy 
effective field theory of QCD, to nuclear force [1], a lot of works have been done during 
last two decades. For reviews, see, e.g., Refs. [2, 3, 4, 5]. Meanwhile, new facilities for 
the rare isotope beams, in which exotic nuclei near the neutron (proton) drip line can be 
created, provide us a new opportunity for investigation [6, 7, 8]. Such a limit close to the 
drip line where the single neutron (proton) separation energy vanishes can be theoretically 
recognized as unitary limit (to be discussed below), and thus the unitary limit would be 
a good theoretical starting point to study the exotic nuclei near the nuclear drip line. In 
this work, we discuss a perturbativc expansion, as the Born series, of an integral equation 
around the unitary limit (as well as weak coupling limit) by sending the cutoff parameter 
to small values. We study it in 5'-wave neutron-deuteron (nd) scattering for spin quartet 
channel in pionless effective field theory (EFT). 

Pionless EFT [9, 10] is one of low energy EFTs for few nucleon systems, in which one 
chooses a typical momentum scale Q smaller than the pion mass, m^, and thus the pions 
are regarded as heavy degrees of freedom and integrated out of effective Lagrangian. Thus 
the large scale of the pionless EFT is A ~ m^r, and the theory provides us a systematic 
expansion scheme in terms of Q/A. In applications of the pionless EFT, much more 
attention has been paid to 5- wave nd scattering for spin doublet {S = 1/2) channel 
because the one-nucleon-cxchangc interaction becomes "singular" [11, 12]. To control the 
singularity one promotes a three-body contact interaction to leading order (LO), and the 
coupling constant of the contact interaction exhibits so called limit-cycle when the scale 
parameter A is sent to the asymptotic limit. This behavior is understood to be associated 
with "Efimov effect". On the other hand, an application of pionless EFT to S'-wave nd 
scattering for spin quartet {S = 3/2) channel is regarded well known because the phase 
shift 6q of the quartet nd scattering is almost perfectly described by two effective range 
parameters in two-body deuteron channel, and no cutoff dependence is reported [13, 14]. 

In this work, we revisit the S'-wave quartet nd scattering below deuteron breakup 
threshold in pionless EFT and reduce the cutoff parameter A smaller than m^. The small 
scale limit might be interesting for studying, e.g., a relation between the pionless EFT 
and a Halo/Cluster EFT [15] in which one may choose a typical scale smaller than the 
deuteron breakup momentum, and the deuteron is never broken into two nucleons and 
could be regarded as a cluster or an elementary field [16]. 

When the value of A is sent to smaller than m^, the cutoff dependence emerges even in 
the quartet channel. Thus we introduce a three-body contact interaction when we solve 
the integral equation using the small cutoff values. After renormalizing the strength of 
the three-body interaction by using the scattering length 04 of the S'-wave quartet nd 
scattering we expand the integral equation for the scattering length and the scattering 
amplitude in terms of the Born series up to next-to-next-to leading order (NNLO) ^. We 
expand the Born series around so called trivial and non-trivial fixed point studied in the 
renormalization group analysis by Birse, McGovern, and Richardson [17]. The trivial 

^We note that this is not an expansion scheme in EFT, but in terms of the Born series. We do not 
have a clear expansion parameter as Q/A. 
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fixed point corresponds to weak coupling limit where all interactions vanish, whereas the 
non-trivial fixed point does to unitary limit where the scattering length becomes infinity 
or the binding energy vanishes. Around the unitary limit, the Born series is expanded 
around the inverse of the amplitude [18]. Wc find that, by reducing the cutoff value, both 
a4 and So are considerably well reproduced by a few terms of the Born series expanded 
around the unitary limit. 

This work is organized as the following. In Sec. 2, effective Lagrangian is displayed. 
In Sec. 3, integral equations for the scattering length and the amplitude are given, and 
the coupling of the three-body interaction is renormalized by the scattering length 04. 
In Sec. 4, the scattering length and the amplitude are expanded in terms of the Born 
series around the weak coupling limit and the unitary limit, and the numerical results are 
obtained. Finally, in Sec. 5, the discussion and conclusions are presented. 

2. Effective Lagrangian 

The effective Lagrangian for the 5'-wave quartet nd scattering in pionless EFT reads [10, 

12] 

C = U + Ct + Cs, (1) 
where jCn is the standard one-nucleon Lagrangian in the heavy-baryon formalism, 

Cn = {iv-D + -^\{vDf-D^]+---]N, (2) 

where is a velocity vector satisfying a condition v"^ — 1, is the covariant derivative, 
and niN is the nucleon mass. Ct is dibaryon effective Lagrangian for ^5"! state. 



a4 ^iv.D+^ [{v ■ Df - D^] + A,} t, - Vt {4 [n^P^'^^n] + h.c] , (3) 



where at is a sign factor, at = ±1. ti is a dibaryon field for spin triplet (^5*1) state. 
is the covariant derivative for the dibaryon field. is the mass difference between the 
dibaryon and two nucleons, = 2mjv + A^. is a coupling constant for dibaryon- 
nucleon-nucleon (dNN) vertex. P/ '^^^ is a projection operator for two nucleons in the ^5"! 
state, 

Pf''^ = ^T^a^a, , (4) 

where Tq and (jj are Pauli matrices for isospin and spin, respectively. Those two constants. 
At and i/t, and the sign factor at {— —1) are determined by two effective range parameters, 
deuteron binding momentum 7 and effective range p^, in NN scattering for "^Si channel. 

We obtain £3, an effective Lagrangian of three-body contact interaction for the 5"- wave 
spin quartet nd state, as 

>C3 = -'^g{A)y^N^[{a-i)\a-^ + 3{a-^{a-t)^]N. (5) 
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Figure 1: Diagrams for 5'-wave nd scattering in spin quartet {S = 3/2) channel. A siiaded 
blob denotes scattering amplitude a line (a double-line) nucleon (dibaryon) , and a double 
line with a filled circle dressed dibaryon propagator (see Fig. 2) . 




Figure 2: Diagrams for a dressed dibaryon propagator. See the caption of Fig. 1 as well. 



We note that the three-body interaction is a higher order term in the quartet nd scattering 
in pionless EFT, and it is not necessary to be introduced when the value of the cutoff 
A is about ttIt^ or larger, A > ttItt, at leading order (LO). We introduce it, however, so 
as to reproduce the scattering length 04 in the quartet channel (and the strength of the 
coupling g{A) is adjusted) when the A value is reduced smaller than uIt^, A < m^r. 

3. Integral equations 

3.1 Integral equation of ^-wave nd scattering amplitude for quartet channel 

Diagrams for iS-wave nd scattering for the spin quartet {S = 3/2) state are given in 
Fig. 1, whereas those for the two-body part in Fig. 1, the dressed dibaryon propagator, 
are given in Fig. 2. The integral equation of the S-wave quartet nd scattering without 
the three-body contact interaction is known as [14] 

t(p,k) = K^''\p,k]E) 

mpd 

[\,K^%,,;E) , ^'^(^''^ , , (6) 

^ -^0 -7 + |Pd (7' - + mE) + g2 - mE 



where t(p, k) is the half-off-shell scattering amplitude and p (k) is the magnitude of off- 
shell final state (on-shell initial state) relative three momentum in CM frame. K^'^^p, q; E) 
is the one-nucleon-exchange interaction. 



T^(a)( 1 , (p^ + +pq-mE\ 

>{p,q;E) = — In ^ -], 7) 

2pq \ -\- — pq — mE 



where E — — ^ + '^[^k'^'- 7 and pa are the deuteron binding momentum and the effective 
range, respectively. A sharp cutoff A is introduced in the loop integral. 
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We note that because a singularity at g ~ 197 MeV which represents an unphysical 
deeply bound state appears in the propagator, the dressed dibaryon propagator is usually 
expanded in terms of the effective range, pd, to avoid an effect from the unphysical bound 
state. However, we are interested in reducing the cutoff value, less than A = 197 McV, so 
we do not expand the dressed dibaryon propagator in the integral equation in Eq. (6). 

Thus the on-shell T-matrix is given by 

T{k,k)^^dt(k,k)^d, (8) 

where is the deuteron wavefunction normalization factor, Z^, = 7Pd/(l — 7Pd)- In terms 
of the half-off-shell scattering amplitude T(p, k), we have the integral equation as 

T{p,k) = -J^ZaK^^\p,k;E) 

rriMPd 



3.2. Integral equation for scattering length 

The scattering length, 04, of the quartet nd scattering is obtained by the on-shell 
scattering amplitude with zero momentum as 

r(0,0) = -— 04. (10) 

TUN 

Now we introduce a half off-shell amplitude (or a half off-shell scattering length) as 

37r 

r(p,0) = a(p,0), (11) 

rriN 

where a(0, 0) = 04. Thus we have an integral equation in terms of a(p, 0) as 
a(p,0;A) = ^X(«)(p,0;-S2) 



--/ dqK^^\p,q;-B,) / / ^ , ^ 3 . a(g,0;A), (12) 



where we have included the A dependence in the scattering length a(p, 0; A) in the above 
expression though a(p, 0; A) is insensitive to A when A ~ ttItt or larger. B2 is the deuteron 
binding energy, B2 = ■j^/ruN- 

Now we choose A = 197 MeV, and solve Eq. (12) numerically. We have two input 
parameters, 7 and pd, whose values are known as 7 = 45.7 MeV and pd — 1.76 fm, and 
thus have 

af = a(0, 0; A = 197 MeV) = 6.34 fm . (13) 
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Figure 3: Left panel: Half off-shell a{p, 0; A) (fm) with A = 197 MeV as a function of off- 
shell momentum p (MeV). Right panel: Cutoff dependence of scattering length a(0, 0; A) 
(fm) without the three-body interaction. 



We reproduce results obtained by Bedaquc, Hammer, and van Kolck [13, 14], 04 = 6.33 ± 
0.10 fm, and Griesshammer [19], 6.354 ± 0.020 fm, up to NNLO in pionless EFT. It agrees 
well with those obtained by Bedaque and Griesshammer [20], 6.8 ± 0.7 fm, up to NLO 
in EFT with perturbative pions, H. Witala et al. [21], 6.321 ... 6.347 fm, from various 
combinations of modern NN potentials and three-body forces, and the experimental 
datum [22], 

af^*- = 6.35 ±0.02 fm. (14) 

In the left panel of Fig. 3, we numerically calculate and plot the half off-shell scattering 
length, a(p, 0;A), at A = 197 MeV as a function of the off-shell momentum p. We find 
that the off-diagonal part of the scattering length quickly decreases and becomes almost 
negligible at p > m,r- This fact may indicate the insensitivity of the amplitude for the 
quartet channel to the value of A when A = or larger, A > TOtt. 

We now reduce the value of A. In the limit where A is sent to zero, we have the on-shell 
scattering length a(0, 0; A) as 

8 1 

lima(0,0;A) = — -~ 19.19 fm. (15) 

A^o ^ ' ' ^ 37(1 -7p<i) ^ ^ 

In the right panel of Fig. 3, we plot the cutoff A dependence of the scattering length 
a(0, 0;A). One can see that there is almost no A dependence at A > 100 MeV, whereas 
the significant cutoff dependence appears in the calculated scattering length a(0, 0; A) 
when the value of the cutoff is reduced less than 7 (= 45.7 MeV). 

3.3 Integral equation with three-body contact interaction 

At the small cutoff region, especially being smaller than 7, the scattering length 
a(0, 0; A) has the significant cutoff dependence. To make the result cutoff independent. 
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Figure 4: Left panel: a(0, 0) (fm) in the cases of three fixed cutoff values A = 27, 7, and 
0.57 as functions of ^'(A) (fm^). A dashed horizontal fine for af' is also included. Right 
panel: g{A) (fm^) which reproduce af- (curve) and three-body bound state with B3 — 
(dashed curve) as functions of A (MeV). 



as discussed before, we introduce the three-body contact interaction g{A) which is sup- 
posed to take account of physics integrated out due to the small cutoff, and thus have the 
integral equation as 



a(p,0) = 



8^ 
3pd 



7^ +p- 

A 







dq [K^''\p,q;-B2)+g{A) 
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r 



a(?,0),(16) 



where we have removed the A dependence from the half off-shell scattering length a(p, 0) 
because the A dependence from the integral should be cancelled by g{h). 

In the limit that the cutoff A becomes large, at A = 197 MeV, the coupling constant 
g{K) should vanish, g{A) 0. On the other hand, in the limit that the cutoff A vanishes, 
one has 



8^ (I 

3 pd \r 



(17) 



where g{Q) = —12.56 ■ ■ • fm^ to reproduce the value of a(0, 0; A = 197 MeV) in Eq. (13). 
Thus we consider a range of g{A), —12.56 < 5'(A) < (fm^), in this work. 

In the left panel of Fig. 4, we numerically calculate from Eq. (16) and plot curves of 
a(0, 0) with three fixed values of A, A = 27, 7, and O.57 as functions of g{A). A dashed 
horizontal line for a^J^ is also included in the figure. One can notice that the curve for 
A = 27 is insensitive to g{A) except for a singular point which corresponds to a three-body 
bound state with zero binding energy, = 0, whereas the curves with A = 7 and O.57 
vary widely and smoothly. 
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In the right panel of Fig. 4, we numerically calculate from Eq. (16) and plot curves of 
g{A) which reproduce af- (curve) and the three-body bound state with — (dashed 
curve) as functions of A. One can see that the three-body bound state with — 
appears when the strength of g{A) becomes stronger than that of g{A) which reproduces 
o^'^ and A is larger than about 60 MeV in the figure. Wc find that the curve of g{A) 
for the three-body bound state with B^ = varies smoothly, whereas that of g{A) which 
reproduces of" has a plateau like shape at the middle of the range of A, A ~ 60 ~ 
100 MeV. Wc use the curve of g{A) which reproduces of* when studying the perturbation 
expansions, the Born series, in the following. 

4. Pert ur bat ive expansions: the Born series 

Now we expand the scattering length and the amplitude in terms of the Born series 
around the weak coupling limit and the unitary limit, as discussed in the introduction. 

4.1 The Born series for the scattering length 

Firstly, we study the expansion in terms of the Born series for the scattering length a^. 
Thus the scattering length a(0, 0) in Eq. (16) is expanded as 



where 



bo 
bi 
b2 



SZ 

04 = a(0, 0) = [6o + + &2 + 
3pd 



-2+9^^ 



I dq 

Jo 

/ dq 

Jo 

X / dq' 
Jo 



^2 _|_ q2 
1 



+ ^(A) 



F{q) 



^2 _|_ g2 

K^^\q,q',-B,) + g{A) F{q') 



(18) 

(19) 
(20) 



72 + q'^ 



+ 9W 



(21) 



with 



3n I 



7 + 



7^ + 



iP.(7 + ^/^ 



7) 



(22) 



This expansion corresponds to that around the weak coupling limit. On the other hand, 
we consider another expansion around the inverse of the scattering length as 



1 



5pd 



i_ 61 _ &2 &| 
bo bl bl + bl 



(23) 



This one corresponds to that around the unitary limit. 

In the left and right panel of Fig. 5, we numerically calculate and plot curves of the 
scattering length 04 expanded around the weak coupling limit and the unitary limit. 
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Figure 5: Left panel: Scattering length 04 expanded around weak coupling limit as func- 
tions of A (MeV). Curves labeled by "LO" are results at leading order, "NLO" up to 
next-to-leading order, and "NNLO" up to next-to-next-to leading order in the both pan- 
els. A horizontal line of af- is also included. Right panel: Scattering length 04 expanded 
around unitary limit as functions of A (MeV). 



respectively, as functions of A. Curves labeled by "LO" are results which include only 
the leading order term, bo, those by "NLO" are results up to next-to- leading order (NLO) 
which include first two terms in the brackets in Eq. (18) and (23), and those by "NNLO" 
are results up to next-to-next-to leading order (NNLO) which include all terms in the 
brackets. 

In the left panel of Fig. 5, one can see that a region of the cutoff A, where the curves of 
obtained from the terms up to NLO and NNLO expanded around the weak coupling limit 
agree with a^J^, is quite small, up to about 10 McV. It would be a natural consequence of 
the perturbation around the weak coupling limit because the perturbation would converge 
when the scale A becomes smaller than a typical scale of the process. In the present case, 
it may be 04 {1/af ~ 31.4 MeV), and thus it converges when the cutoff A becomes much 
smaller than 1/04'*. In the right panel of Fig. 5, on the other hand, we find that the region 
where a*/*^ is reproduced is broaden, up to about 100 McV due to the two or three terms 
(up to NLO and NNLO, respectively) of the Born series expanded around the unitary 
limit. 

4.2 The Born series for the scattering amplitude 

Before expanding the integral equation for the scattering amplitude in terms of the 
Born series, we check how the three-body interaction g{A) can reproduce the phase shift 
Sq at the small cutoff values. In Fig. 6, we numerically calculate and plot the phase shift 
Sq below deuteron breakup threshold as functions of the on-shcll momentum k by solving 
the integral equation at small cutoff values, A = 197, 140, 100, and 60 McV. The deuteron 
breakup momentum is kbr — 52.7 MeV. We also include results from a modern potential 
model (Avl8) [23] in the figure. One can see that the results are fairly independent on 
the cutoff values and considerably well agree to those obtained from the accurate Avl8 
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Figure 6: Phase shift Sq (deg.) of the S'-wave nd scattering for spin quartet channel below 
deuteron breakup threshold as functions of k (MeV). Curves are obtained by solving 
the integral equation at cutoff values, A = 197, 140, 100, and 60 MeV with the three- 
body interaction g{A) renormalized by af-. Plus signs "+" denote results from a modern 
potential model (Avl8) [23]. The deuteron breakup momentum is kbr — 52.7 MeV. 



potential model. 

We now expand the on-shell T-matrix in terms of the Born series up to NNLO around 
the weak coupling limit and the unitary limit. Thus the on-shell T-matrix is obtained as 

8nZd 



T{k,k) = 



rriNPd 



[Bo + B1 + B2 + 



with 



where 



Bo = 
Bi = 

Bo = 



Giq,k) 



V{k,k)^K^^\k,k-E)+g{K), 

[''dqV{k,q)G{q,k)V{q,k), 
Jo 

r dqV(k, q)G(q, k) T dq' V{q, q')G{q\ k)V(q', k) , 
Jq Jq 



8 



7 + \h' + iii'-k') 



(24) 

(25) 
(26) 

(27) 
(28) 



1 - Ip, + ^72 + |(g2_p)J g2 _ fc2 _ 

This corresponds to the expansion around the weak coupling limit. We also have an 
expansion around the inverse of the T-matrix as 



1 



T(k,k) 



1 Bi B2 B^ 

Wo'bI'bI^bI^ 



(29) 
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This expansion corresponds to that around the unitary hmit. 

To calculate the phase shift Sq from the Born series around the weak coupling limit in 
Eq. (24) we employ the relation 



l + z^—T(k,k) 



(30) 



We note that when the ampHtude expanded in the Born series is truncated, it breaks 
unitary condition and the phase shift becomes a complex number [24]. On the other 
hand, to calculate the phase shift from the expansion around the unitary limit in Eq. (29) 
we use the formula 



kcot Sq — ^^Re 



ruN T{k,k)' 



with 



Re 
Re 



1 



TNLo{k,k) 
1 



TNNLo{k, k) 



1 

1 



ReBi 

ReBi 
Bn 



ReB2 , (i?e5i)2 



5n 



+ 



{ImBif 



Bl 



(31) 

(32) 
(33) 



where Bq has real part only. This expansion preserves the unitary condition, at least up 
to NNLO. 

In Fig. 7, we numerically calculate and plot real and imaginary part of the phase shift 
5q of 5'-wave quartet nd scattering obtained from the truncated Born series expanded 
around the weak couphng hmit in Eq. (24). Curves up to NLO are obtained by including 
first two terms, Bq and and those up to NNLO by including all three terms, -Bq) Bi, 
and i?2, in Eq. (24). We have fixed A at A = 13 MeV because we found in the left panel 
of Fig. 5 that the scattering length a}^ are fairly well reproduced when the cutoff value 
is about up to 10 MeV. A line labeled by "FuU"^ obtained from the calculation without 
the truncation at A = 197 MeV is also included in the figure. One can see that the 
real parts of Sq agree with the full result up to about 6 McV for NLO and 9 McV for 
NNLO. In addition, the imaginary parts sharply decrease around the edge of the cutoff 
value A = 13 MeV, and it indicates that the unitary condition is broken. The broken 
unitary condition up to NLO is also partly cured by including the higher order term, i?2, 
at NNLO. 

In Figs. 8 and 9, we numerically calculate and plot curves of the phase shift 5q obtained 
from the Born series up to NLO and NNLO in Eqs. (32) and (33), respectively, expanded 
around the unitary limit where the cutoff values are chosen A = 140, 100, and 60 MeV. 
In addition, plus signs "+" labeled by "Full" in the figures denote the results from the 
calculation without the truncation at A = 197 MeV. We find in Fig. 8 that the curves 
considerably well converge to that of the full calculation as the cutoff value decreases. 



^This line is the same as that at A = 197 MeV in Fig. 6. 
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Figure 7: Real and imaginary part of phase shift 6o (deg.) of S'-wave quartet nd scattering, 
as functions of k (MeV), obtained from the truncated Born series up to NLO and NNLO 
expanded around the weak couphng hmit where the cutoff is fixed at A = 13 MeV. A 
dashed hne labeled by "Full" obtained from the calculation without the truncation and 
A = 197 MeV is also included. 



and the truncated Born series expanded around the unitary limit up to NLO fairly well 
reproduces the result of the full calculation at A = 60 MeV. In Fig. 9, we can see that 
the convergence to the full result, as reducing the cutoff value, becomes faster due to the 
inclusion of the higher order terms. 

5. Discussion and conclusions 

In this work, we studied the perturbative expansions, as the Born series, of the integral 
equation around the weak coupling limit and the unitary limit for the S'-wave nd scattering 
for the spin quartet channel below the deuteron breakup threshold in pionless EFT at the 
small cutoff values. The three-body contact interaction is introduced when the integral 
equation is solved by using the small cutoff values. After the strength of the three-body 
interaction is renormalized by using the scattering length 04, we expand the integral 
equation for the scattering length and the amplitude, as the Born series, up to NNLO 
around the weak coupling limit and the unitary limit. We find that the scattering length 
and the phase shift are considerably well reproduced by a few terms of the Born series as 
we reduce the cutoff values to about 10 MeV for the expansion abound the weak couphng 
limit and to about 60 MeV for that around the unitary limit. 

In the present particular process, the S'-wave quartet nd scattering studying, one may 
regard that there is no advantage by sending A to small values because the physical 
observables, the scattering length 04 and the phase shift Sq, are well described (without 
fitting any unknown parameters) by the two effective range parameters in the deuteron 
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Figure 8: Phase shift Sq (deg.) of S'-wave quartet nd scattering, as functions of k (MeV), 
obtained from the Born series up to NLO expanded around the unitary hmit with A = 140, 
100, and 60 MeV. Plus signs "+" labeled by "Full" are obtained from the calculation 
without the truncation at A = 197 MeV. 



channel at the usual large scale of the pionless theory, A ~ m^r. However, it could be a 
useful limit when one studies a relation between pionless EFT and a Halo/Cluster EFT 
whose large scale smaller than m^, such as a deuteron cluster theory for the nd scattering 
below the deuteron breakup momentum. 

In addition, it may be an interesting observation that the physical observables can be 
well described by a few terms of the Born series expanded around the unitary limit at the 
small cutoff values. If this property were common in some class of the interactions and/or 
in the unitary limit, it could provide us a useful method to make a non-perturbative 
interaction perturbative and/or to be used in studies of the exotic nuclei near the drip 
line. Indeed, one may regard that what we found in this work can be a coincidence, a 
unique property only in the S'-wave quartet nd scattering or in the case of "non-singular" 
interactions. Therefore, it may be worth studying whether the same property which we 
found in this work can be observed or not in the case of the singular interactions, such as 
the S'-wave spin doublet nd scattering or the NN scattering involving the tensor force in 
the ^Si-^Di channel. 
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Figure 9: Phase shift 6q (deg.) obtained from the Born series up to NNLO expanded 
around the unitary hmit with A = 140, 100, and 60 MeV. See the caption of Fig. 8 as 
welL 
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